Binary black holes' effects on electromagnetic fields 
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In addition to producing gravitational waves (GW), the dynamics of a binary black hole system 
could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and 
electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected 
by a pair of orbiting black holes through the merger. In particular, we show how the binary's dy- 
namics induce a variability in possible electromagnetically induced emissions as well as a possible 
enhancement of electromagnetic fields during the late-merge and merger epochs. These time depen- 
dent features will likely leave their imprint in processes generating detectable emissions and can be 
exploited in the detection of electromagnetic counterparts of gravitational waves. 



Introduction: As the era approaches when the detec- 
tion and analysis of gravitational waves joins electromag- 
netic observations, studying a number of astrophysical 
systems in both wavebands becomes an exciting prospect. 
Indeed, most astrophysical systems which produce grav- 
itational waves are very likely also electromagnetically 
bright (see e.g. QUI]), with the exception of systems com- 
posed solely of black holes. 

One example of such a system is a black hole sur- 
rounded by an accretion disk. This scenario has been 
used to explain a vast range of spectacular phenomena, 
from AGNs and quasars to gamma ray bursts, all of which 
are prodigious emitters of electromagnetic energy. Strong 
emission from these systems is understood as the result 
of radiative processes within jets powered by the extrac- 
tion of rotational and binding energy. 

The pioneering work of Penrose [3[ and Blandford and 
Znajeck [4|, together with a large body of subsequent 
work, has provided a basic understanding of possible 
mechanisms to explain highly energetic emissions from 
single black hole systems interacting with surrounding 
plasmas. The interaction of electromagnetic field lines 
with the strong gravitational field of a rotating black hole 
is the fundamental component of these mechanisms to ex- 
plain the acceleration of particles that traverse the black 
hole's ergosphere. The scenario of a pseudo-stationary, 
single black hole interacting with an accretion disk is rea- 
sonably well understood, and this system is employed to 
explain energetic phenomena such as gamma ray bursts, 
AGNs, quasars, blazars, etc. However, a highly dynam- 
ical stage may occur prior to such a pseudo-stationary 
regime which could give rise to strong emissions. 

A related system, studied here, involves the collision of 
super massive black holes, whose gravitational radiation 
would be detectable by the Laser Interferometric Space 
Antenna (LISA). To date, strong observational evidence 
indicates that massive black holes exist in the centers of 



most galaxies, and that galaxies have undergone mergers 
in their past. In galactic collisions, the initial black holes 
eventually merge as their orbits shrink through a variety 
of mechanisms. As discussed in [H|, 0, a circumbinary 
disk is formed as the binary hollows out the surround- 
ing gas (see, e.g. [1, 0, [Toj ) • As the distance between the 
black holes decreases, the dynamics of the binary become 
governed by the emission of gravitational waves and dis- 
connected from the properties of the disk [13] • 

While the late time dynamics of the system can be un- 
derstood in terms of the pseudo-stationary picture men- 
tioned above, an interesting intermediate regime would 
include magnetic fields, anchored at the circumbinary 
disk, being influenced by the orbiting behavior of the 
two black holes. In particular, both precursor emission 
and enhancement of electromagnetic energy could be pos- 
sible as the electromagnetic fields interact with the gas 
being stirred and compressed by the black holes' shrink- 
ing orbit. Consequently, the late pre-merger and merger 
dynamics of a supermassive binary black hole system 
could have a strong influence on the possible interac- 
tion between the remaining gas and the electromagnetic 
fields and even seed several emission mechanisms more 
strongly. In this work we study such a scenario by con- 
sidering the Einstein-Maxwell system in a setup that de- 
scribes a pair of black holes close to the merger epoch, 
and we examine the electromagnetic field behavior and 
compare to the single black hole case. While we do not 
consider a plasma in our current study, our analysis helps 
to understand the possible behavior in its presence and 
lay the foundations to future work in this direction. 
Overview of the numerical approach: We solve the 
coupled Einstein-Maxwell system to model the black hole 
merger interacting with an externally sourced magnetic 
field. The Einstein equations are written in the General- 
ized Harmonic formalism described in [13]. The Maxwell 
equations are written directly in terms of the electric 



2 



and magnetic fields, as in [14|. The electromagnetic 
fields have astrophysically motivated magnitudes, and 
thus their energy is several orders of magnitude smaller 
than the gravitational field energy. Hence the dynami- 
cal effects of the electromagnetic fields on the geometry 
are negligible. This fact is corroborated by simulations 
not displaying any difference in the black holes dynam- 
ics with an initial magnetic field an order of magnitude 
larger and even zero, than the one employed here. 

In both systems of equations, the constraints are 
kept under control via damping mechanisms: constraint 
damping for the Einstein equations an d divergence 
cleaning for the Maxwell equations (similar to that de- 
fined in [l6[ for ideal MHD). A combination of Sommer- 
feld and constraint preserving boundary conditions are 
applied at the outer boundary [17] for both systems. The 
incoming modes of the electromagnetic fields are defined 
via Dirichlet (maximally dissipative) conditions induced 
from the physical picture of a circumbinary disk present 
beyond the computational domain. 

We use the HAD computational infrastructure that pro- 
vides distributed Berger- Oliger style Adaptive Mesh Re- 
finement (AMR) [11, El with full sub-cycling in time, to- 
gether with a novel treatment of artificial boundaries [20[ . 
The refinement regions are determined using truncation 
error estimation via a shadow hierarchy. A fourth or- 
der spatial discretization satisfying a summation by parts 
rule together with a third order Runge-Kutta scheme for 
the time integration are used to help ensure stability of 
the numerical implementation [2lj |. 

Overview of the physical set-up: To explore the ef- 
fects of the merger dynamics on the electromagnetic field, 
we compare certain cases of a single spinning black hole 
with cases of merging holes. In all cases, the orbital plane 
of the evolution (or equatorial plane for the single BH) 
is assumed to be aligned with that of the circumbinary 
disk. The magnetic field is defined as anchored in the 
disk; hence, its associated magnetic dipole is aligned with 
the orbital and spin angular momentum. 

As we are primarily interested in the dynamics close 
to the merger, we adopt initial data such that it takes 
place after about one orbit. These data describe quasi- 
equilibrium, equal-mass, non-spinning black holes con- 
structed by the publicly available lorene code [22]. The 
black holes have masses given by M s = M/2, and are ini- 
tially separated by « 6M, lying beyond the approximate 
inner most stable circular orbit (ISCO) [23[. The ini- 
tial magnetic field is poloidal and constructed from the 
electromagnetic potential produced by a circular loop, 
whose radius is assumed to be larger than the domain 
of interest Q. We assume the disk lies at ~ 10 3 M 
and so B l ~ B Q z. The electric field is initially zero 
throughout and the magnetic field strength adopted is 
B = 1O 4 (M/1O 8 M ) G, which is within values inferred 
in relevant astrophysical systems (see e.g. [Hj). 

We adopt a cubical domain given by [—125 M, 125 M] 3 



and employ an AMR configuration with 7 levels of refine- 
ment that adjust themselves dynamically to ensure that 
the solution's error is below a pre-determined threshold 
using a shadow hierarchy. Residuals of the various con- 
straints and convergence of the physical fields with reso- 
lution are checked as tests of the code. 
Results: For concreteness we present the results of a sin- 
gle black hole with spin parameter a = 0.7 /M s . We adopt 
this value, which is slightly larger than the final spin ex- 
pected for a merged black hole [34| for comparison pur- 
poses. The evolution shows an initial transient, where the 
EM fields adapt to the geometry of the black hole space- 
time, giving rise to an electric field and a deformation of 
the magnetic field (see also jifl Hi}). After t ~ 40 - 60M 
the solution is clearly seen to evolve towards a quasi- 
stationary state determined by Wald's solution [28[ for 
a Kerr black hole immersed in a uniform magnetic field 
which is aligned with its spin. Fig. [T] presents the electric 
field obtained at t = 200M in the plane y = 0. For com- 
parison purposes, the corresponding field from Wald's 
exact solution is shown for x > 0. The apparent agree- 
ment along with a careful examination of the asymptotic 
solution indicates that, for a black hole immersed within 
an almost uniform magnetic field aligned with its spin, 
the final state is Wald's solution 12611 . 
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FIG. 1: The electric field for y = 0, z > at t = 200M for a 
single black hole, together with the apparent horizon (green) 
and the ergosphere (magenta). The region x < display the 
numerical solution, while Wald's exact solution in x > 0. 

Next, we consider a binary black hole system at a 
quasi-circular stage representing a late inspiral. In the 
present work, we adopt initially non-spinning equal mass 
black holes to explore the dynamics of the system in a 
simple configuration. The initial electromagnetic field is 
defined in exactly the same manner as in the single black 
hole case. The evolution of this system, however, is signif- 
icantly different as the orbiting black holes stir the fields 
during the inspiral. Examples of the field configurations 
during the evolution are provided in Fig. [2j 

To analyze the influence of the binary's dynamics on 
the electromagnetic field we monitor two particular quan- 
tities. One is the electromagnetic energy density E 2 +B 2 , 
while the other is the radial component of the Poynt- 
ing vector, S r (radial from the origin). The former is 
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employed to get a sense of the energy variation in the 
EM field, while the latter is used to examine the induced 
multipolar structure and for comparisons to the gravita- 
tional wave signal computed from the Newman-Penrose 
scalar, ^4. (We additionally compute the analogous EM 
Newman-Penrose radiative quantity $2 and confirm that 
the same features are present in it as in the Poynting 
vector). Note that since the magnetic field is anchored 
in the disk, the magnetic field does not decay with the 
distance from the binary, and this obscures the standard 
interpretations drawn from the Poynting field. Neverthe- 
less, the field does serve to elucidate the mode structure 
induced by the dynamics. 

An examination of these quantities reveals several dis- 
tinct features during the evolution. After an initial tran- 
sient signal, and well before the merger takes place, the 
electromagnetic fields display a pattern consistent with 
those produced by electric dipoles, pointing along the or- 
bital plane in the direction orthogonal to their velocity, 
orbiting about each other. This can be understood from 
the 'membrane paradigm' point of view [29[, in which 
the stretched horizon is endowed with a surface density 
of electrical charge and a surface resistivity. The quasi- 
circular trajectories of the black holes cause a charge sep- 
aration in the direction perpendicular to both the veloc- 
ity and magnetic field, as in the Hall effect. In particular, 
the evolution induces a toroidal electric field (E^) with 
strength given by i^orbitaic -1 . These fields would seed 
a Blandford-Znajek type mechanism and generate energy 
when interacting with a surrounding medium. Addition- 
ally, the electric and magnetic fields oscillate with a fre- 
quency equal to the dominant frequency of the gravita- 
tional radiation. Since the Poynting flux is determined 
by the square of the electromagnetic fields, its oscillation 
pattern has a frequency that is twice that of the domi- 
nant gravitational wave signal, and can be described in 
terms of / = 0, rn = and I = 2, m = modes. As the 
merger approaches, an I = 2,m = ±2 mode is induced 
and the peak of EM radiated energy lags behind the grav- 
itational one by t ~ 2M (1 x 1O" 5 (M/M ) s ) with the 
same frequency (see Fig. [3]). This lag time is likely re- 
lated to the transition time from the binary black hole 
system to a single black hole, measured as the 'time of 
flight' from the ISCO to the light ring. Indeed estimates 
of the ISCO and the light ring (lr) for a = 0.67 jM reveals 
A = nsco — Hr — 1.5M. Additionally, an enhancement 
of the electromagnetic field energy is clearly visible in 
the binary case, which displays a ~ 30 — 40% increase 
with respect to the single case, as illustrated in Fig. [H 
(We conjecture that in the presence of gas and black hole 
spins, this number could significanly increase). At late 
times, the system asymptotes to the same solution ob- 
served for the single black hole. 

Final Comments: We have analyzed the behavior of 
electromagnetic fields influenced by the dynamics of a 
binary black hole system. Our study illustrates several 




FIG. 2: Magnetic and electric field lines at t ~ -40M, -20M 
in 3D (the merger happens at t = 0). The electric field lines 
are twisted around the black hole, while the magnetic lines 
are mostly aligned with the z-axis. 
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FIG. 3: The radial component of the Poynting vector for 
the binary black hole case. The structure reveals it is dom- 
inated by (spherical harmonic) EM co,o and EM C2,o modes in 
the orbiting stage though close to merger a EM 02,2 is induced 
which lags behind the corresponding gravitational wave (spin- 
2 spherical harmonic) GW C2,2- The transient behavior before 
t ~ —100 is spurious and due to the initial data employed. 



interesting aspects of such systems that emit not only 
gravitational waves, but can also radiate electromagnet- 
ically when interacting with a plasma. 

The EM fields have a clearly discernible pattern tied 
to the dynamics, making them additional tracers of the 
spacetime, as these features would imprint particular 
characteristics in processes producing observable EM sig- 
nals. In particular, in the pre-merger stages, the black 
hole dynamics induce EM energy flux oscillations with 
a period half that of the dominant GW signal produced 
by the system, i.e., a fourth of the orbital period, and 
a gradual enhancement of the energy in the electromag- 
netic field. During the merger, additional structure is 
impressed on the field, yielding a distinct pattern that 
lags behind the gravitational wave output by a time con- 
sistent with the characteristic time from the ISCO to the 
light ring of the merged black hole [35[ . 

The distinct variability of the electromagnetic field 
as it is dragged by the black holes suggests tantalizing 
prospects to detect pre-merger electromagnetic signals 
from systems detectable in the gravitational wave band. 
At a more speculative level, these combined signals could 




FIG. 4: Normalized electromagnetic energy density (with re- 
spect to its value at the initial time in the absence of black 
holes). Left/Right columns corresponds to y = and z = 
plane while the first three rows correspond to times t ~ —70, 
5 and 30M and the last to the single black hole case at 
t — 150M for the sake of comparison. Normalized energy 
contour plots are shown with lines from 1.2 to 0.8 at intervals 
of 0.1. Around the merger epoch, the energy computed for 
the binary case is ~ 30 — 40% larger than in the single case. 



be exploited to study alternative theories where photons 
and gravitons might propagate at different sjjeeds (for a 
recent discussion of some possibilities see [30|, El). A 
complete description of the problem requires the incor- 
poration of gas and radiation effects. However, the main 
qualitative features — driven by the orbiting behavior of 
the black holes, whose inertia is many orders of magni- 
tude above all else — would intuitively remain unaltered. 
The same reasoning also indicates the qualitative fea- 
tures described above will be present in non-equal mass 
binaries. However, the presence of spin could introduce 



further features, as energy could be extracted from the 
individual black holes prior to merger. Ongoing [32[ and 
related work will help to test these conjectures and shed 
further light on the problem. 
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